Fractional exclusion statistics — the method to describe interacting particle systems as 

ideal gases 
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I show that if the total energy of a system of interacting particles may be written as a sum of 
quasiparticle energies, then the system of quasiparticles can be viewed in general as an ideal gas 
with fractional exclusion statistics (FES). The general method for calculating the FES parameters 
is also provided. 

The interacting particle system cannot be described as an ideal gas of Bose and Fermi quasipar- 
ticles except in trivial situations. 
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INTRODUCTION 

A fractional exclusion statistics (FES) system [U-El 
consists of a countable number of species. I will count 
the species using the indexes i and j. Each species con- 
tains Gi single-particle states and Ni particles. The FES 
character of the systems consists in the dependence of the 
GVs on the A,'s. For example by adding SNi particles 
in the species i the dimensions of all the other species 



change according to 6Gj 



-ctjiSNi, for any i and j. 



The parameters are called the FES parameters. 

In practice, the species may be different types of par- 
ticles that coexist in the same system, quasiparticle exci- 
tations in the lowest Landau level in the fractional quan- 
tum Hall effect EL J2 , excitations [l| or motifs of spins 
in spin chains [J, |5j, elementary volumes obtained by 
coarse-graining in the phase-space of a system [3, 
etc. Some amendments to the FES formalism have been 
introduced more recently, in Refs. @,@"11|. 

In this paper I shall use the ansatz 

ctij = G l a lJ + ai6ij. (1) 

which applies quite generally in quasi-continuous systems 
1- 

The equilibrium distribution of particles in the species 
may be calculated in two equivalent formulations: I call 
them the Bose and Fermi perspectives Q . 

In the Bose perspective I use Gi to denote the number 
of available single-particle states of species i. In such a 
case, the number of microconfigurations in which we can 
find the system is W = Ui{( G i + N i ~ 1) ( G i ~ 1) ! ] }• 
If we add a small perturbation, {SNi}, to the particle 
distribution, W becomes 



W B = \ 



[d + Ni — 1 + (1 — a z )SNi - £V OjjSNjY. 
{N t + SNi)\(Gi - 1 - a t SN q - £. a l0 5N 3 )\ ' 



(2) 

Associating the energies {e^} to the particles in each of 
the species, we write the partition function, Z 0], which 



we then maximize to obtain a set of equations for the 
equilibrium particle distribution, 



00* 



In- 



1 ' • ( 3 ) 



where (3 = l/(fcflT), \i is the chemical potential of the 
system, and hi = Ni/Gi is the equilibrium particle pop- 
ulation. 

In the Fermi perspective, T, = Gi + Ni — 1, and 
Wf = Udn/WiTi - Aj)!]}. At the variations {SNi}, 
we obtain 



W F 



n 



{Ti-aiSNi-Gi^a^SNj)] 



[Ti - ^ - (1 + ai )SNi - Ti £ . oyflV 3 -]! 
1 



(4) 



(Ni + 5Ni)\ 

and the equations for the equilibrium particle popula 
tions, ft = Ni/T h are 

[1 - h] 1+a 



P(m - e») + In 



fi 



X'.-"- 1 " 1 ( 5 ) 



I will show below that systems of interacting particles 
may be described as ideal FES systems. For interacting 
bosons, more natural is the Bose perspective, whereas for 
interacting fermions it is more convenient to employ the 
Fermi perspective. 

Equations ([3]) and ([5]) can be readily transformed into 
integral equations in the quasicontinuous case. If in- 
stead of the index i (or j) we introduce the quasicon- 
tinuous variable i (or j) - which may be a multidi- 
mensional variable, like the quasi-momentum, or a one- 
dimensional variable, like the quasi-energy - of density 
of states (DOS) <j(i), then Eqs. © and © become [U 

/3( Ati -e i )+ln [1 + y 1 - = /V(j)ln[l + & j ]a ji dj. (6) 



and 
P(fH-ei) 



-In 



[i-Zt] 1 -" 1 

fi 



cr(j)ln[l-/j]ajidj, (7) 
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respectively. 

THE QUASIPARTICLES 

The ideal gas description 

Suppose we have a system of interacting particles 
which we want to describe as an ideal gas of quasipar- 
ticles. For this we introduce the quasiparticle energies, 
{fii}, which we want to satisfy three conditions, specific 
to ideal gases: 

• E({m}) = ^riiii, (8a) 

i 

• the energies {ii} are well defined and therefore 
independent of the set of occupation numbers, 

{m}; (8b) 

• the equilibrium populations, (rij)(T, fi, e,), are 
functions of only T, /x and ii and are independent 
of the populations. (8c) 

The heat capacity of any system (in units of k B ) is 



C v fdU 



k B \8TJ N \8T /fi 



dU\ ( dU\ fdN\ fdN 



d^J T \dT / 



If the conditions ([SJ arc satisfied, then 



dU 



d{k B T) 



(SIP 
\dfj, /7 



dN 



dfj, ) T 
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E 



d(k B T)^ 
d( ni ) \ 

( djm) \ 
\d(k B T))> 

d(nj) 
d/x 



(9) 



(10a) 
(10b) 
(10c) 
(lOd) 



where we used the notation U(T,n) = (E)x^ for the in- 
ternal energy of the system, at temperature T and chem- 
ical potential, /x. 

For an ideal gas of fermions in the quasicontinuous 
limit and with a density of states (DOS) of the form 
a(i) = Ce s , with C and s being two constants 12, 13j . 
the total particle number and the internal energy are 



N 



C{k B T) s+L T(s + l)Li s+1 



and (11a) 



U = C(k B T) s+2 T( S + 2)Li s+2 (-e^), (lib) 

respectively, where the function Li n is the polylogarith- 
mic function of order n 1J, [15 1 . 



Using Eqs. (|T0|) we calculate the specific heat (also in 
units of k B ): 



k B 



Cy 
k R N 



(« + !)(« + 2) 



Li 



s+2 



Li s+1 (-e^) 



-(■s + l) 2 



s+l 



(11c) 



In the low temperature limit /3/i 3> 1 and in the lowest 
orders of approximation, Eqs. (TTTj) become 



N = C f i s+1 /(s 

U = C[i s+2 /(s + 2) 

Cy 7T 2 k B T 

— w (s + l)- 

k B 6 fj. 



1) ee C~e s F +1 /(s + 1), 



(s + l) : 



(12a) 
(12b) 

(12c) 



where Eq. (|12a[) defines the Fermi energy, ip = [(s + 
l)N/C] x /( s+1 \ and Eq. §FM may be put into the low 
temperature universal form jla . 1 1 7 1 , 



Cy 
k B 



Ncy 7T 2 

-, — = —k B Ta{e F ). 
k B S 



(13) 



The quasiparticle gas 

The gases of quasiparticles used to describe systems of 
interacting particles do not satisfy in general the condi- 
tions ©. An example shown in Ref. 18| is the Lan- 
dau's Fermi liquid theory (FLT). The quasiparticle en- 
ergy, ej, depends on the occupation of the other quasi- 
particle states, i.e. = (k({nj}), where {rij} denotes 
the set of all occupation numbers. In such a case the 
condition ()8bj) is not satisfied and through the pop- 
ulation (nj)(T, /z, ?i) depends on the populations of the 
other quasiparticle levels, violating the condition (f5cj) as 
well. Moreover, in the FLT the sum of the energies of 
the quasiparticles is not equal to the total energy of the 
system, violating also condition (|5aj) (see [H|). 

Let us now see how the heat capacity of the system 
can be calculated. For this I will suppose that condition 
(|8a|) is true, otherwise the gas of quasiparticles may not 
be used for this purpose [18j. 

If ii is a function of {rij}, then Eqs. (|10a[) and (jlObj) 
are not valid because ii varies with T and /x: 



dej(T,(i) 
d{k B T) 
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dej(T,ri d(n )(T,ii) 
d{n 3 ) 8(k B T) [ 

d~e i {T,ii)d{n j ){T, l j) 



d{nj) 



(14a) 



(14b) 



On the other hand, the derivatives of the populations 
are given by the equations 

d{rii){T,n) d(ni)(T,fi,ii) d(rn)(T, (i, e,) dei(T, (j.) 



d(k B T) 



d(k B T) 



d(k B T) ' 
(15a) 



3 



d{rii){T,ij) d{n t )(T 7 p,h) 9(n»)(T, p 7 h) dh{T, p) 



dp 



dp 



dh 



dp 
(15b) 

Plugging Eqs. (| 14|) into (|T5|) we obtain a self-consistent 
system of equations for the variation of the occupation 
numbers with T and p. 

Using the solutions of Eqs. (fT5j) together with Eqs. 
(fig)), we can calculate dU/d(k B T), dU/dp, dN/d(k B T), 
and dN/dp, to finally obtain the heat capacity ([9]), where 
Eqs. (|10aj) and (|10b|) become 



dU 



d{k B T) 



E 



dU 
dp 



E 



dp 



d{k B T) yi 

( dh 
\d(k B T) 



dp / T 



(16a) 



(16b) 



THE IDEAL FES GAS 

Assuming that we choose the quasiparticle energies, h, 
in such a way that condition (|8a|) is satisfied, let us see 
how we can satisfy also the conditions (|8b[) and (|8cl) . 



The principle of the method is given in Refs. 0MIH- 

I change to the quasicontinuous description, i — >• i, and 
I assume for simplicity that i is a ID quantity (e.g. the 
energy of the free particles, e 0,13 Hi)- The quasiparti- 
cle energy, h({n>j}), is then a functional of the occupation 
numbers, {nj}, or, for the equilibrium distribution, it is 
a functional of the populations - 

I choose i or h in such a way that if i < j, then 
h({(nk)}) < ^({("k)})- hi this way I establish a bi- 
jectivc correspondence between i and e, which one may 
invert and write U{{(nj)}) or ig({(ng/)}). 

If <r(i) is the DOS in the variable i, then 



a(e) = «r(i) 



di 




de 








de 


di 



(17) 



a(e) is a functional of With the aid of cr(i) and 

cr(e) we define the particle densities, p(i) = a(i)(n$) and 
p(e) = a(e){m). 

Now I transform the quasiparticle gas into an ideal gas 
by simply changing the perspective: the usual perspec- 
tive is to see the quantum numbers of the ideal gas, i, as 
fixed and the quasiparticle energies, h({(nj)}), as func- 
tionals of the populations (Fig. [T]a). When we invert 
the perspective we see the quasiparticle energies, e, fixed 
and the free particle quantum numbers as functionals of 
the populations, ig({(ng)}) = i?[p(e')] (Fig. d]b). 



Normal perspective: 

i and cr(i) are fixed 
£j[P(i)]is a functional 

Bose or Fermi 
statistics 
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FIG. 1: (Color online) The two perspectives: (a) the the free-particle quantum numbers, i, are fixed and the quasiparticle 
energies ?i[p(j)] are functionals of the particle density, p(j); (b) the quasiparticle energies, e, are fixes and the free-particle 
quantum numbers, ii[p(e)], are functionals of the particle density, /5(e). 



By the relation (JTTj) and since cr(i) is a quantity that is 
fixed by the properties of the single-particle states, cr(e) 
becomes a functional of the particle density, p{e). This 
property ensues FES. 

To show this I coarse-grain both axes, i and e. Each 



interval, [i m ,i m +i] on the i axis or [fm, fm+i] on the 
e axis, represents a species, with G m = c(i m )(i TO+ i — 
im) = cr(e TO )(e TO+ i - hn) and N m = p(i m )(i m +i - i m ) = 
p(e m )(hn+i — hn) (assuming that the intervals are small 
enough to use the linear approximation); I use the letters 



4 



m and n to designate the species. 

In the FES perspective, the insertion of 8N m particles 
in the species m causes a change of the interval [i„, i n +i] 
on the i axis by 



5N„ 



Si 



n+l 



8p{e m ) Sp( 



SN m (i 



m+1" 



l) di 



Sir, 



Sp{e r , 



(18) 

where by Si n /8p(l m ) I denote the functional derivative, 
which I assume to be analytic. 

The change of the interval [i„, i n +i] leads to the change 
of the number of states in the species n by 



SG„ 



SN„ 



— — er(i„+i) - s (r{hi) 
op(e m ) op(e m ) 



SN m (i m+ i - — 
ai 



Sin 

5p(e m 



-o-(i n ) 



(19) 



Mapping back, i onto e, I can express the change of the 
number of states in the species n as 



SG„ = 5N m (e n+ i-e n ) — 



Sp(e m ) 



- a s„e m SN m , 
(20) 

where by 5e n /5p(e m ) I denote the functional derivative 
of e„ with respect to the variation of 5p(e m ), where i is 
fixed - this calculation trick does not change the nonin- 
teracting character of the quasiparticles. 

The last part of Eq. (|20| gives us the FES parameter, 
ag„g m . We notice that az n e m oc (e n +i — e n ) = G n /a(e), 
so it is proportional to the dimension of the species that 
it acts upon, in accordance with the ansatz ([T]) 0, @| and 
with the general FES rules of Rcf . Q . 

If the function [Se n /8p(e m )]d , (e n ) is singular in some 
point, e„, then c%e n e m is n ot proportional to the dimension 
of the species G n and 



SG„ 



5N„ 



Sin 



+ 1 



Sp(e m ) 



cr(i„+i) 



Sir, 



Sp{e m ) 



<y{i n ) 



The parameters a^ s \ obey the rules of Ref. [9[ and in 
addition might still obey the ansatz (JTJ) proposed in the 
Introduction. Such a case was discussed in Ref. @, EH , 
in relation to the Fermi liquid theory. 

Once the FES parameters are known, the thermody- 
namics follows according to the formalism outlined in the 
Introduction. 



(21) 



CONCLUSIONS 

I showed that if a system of interacting particles is 
described as a gas of quasiparticles of energies, e({nj}), 



such that E({m}) — (condition [8a)l . then the 

gas of quasiparticles may be viewed as an ideal gas that 
obeys fractional exclusion statistics (FES). If the condi- 
tion (|8ap is not satisfied, as it happens with the Landau's 
quasiparticles in the in the Fermi liquid theory, then the 
quasiparticle gas cannot be used for the calculation of 
the thermodynamic properties of the original interacting 
particle gas, except in some trivial cases (see fl8| ) . 

The general method for calculating the FES parame- 
ters of the quasiparticle gas is also given. 
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